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We investigate Hawking evaporation in a recently suggested picture in which black holes are Bose 
condensates of gravitons at a quantum critical point. There, evaporation of a black hole is due to two 
intertwined effects. Coherent excitation of a tachyonic breathing mode is responsible for the collapse 
of the condensate, while incoherent scattering of gravitons leads to Hawking radiation. To explore 
this, we consider a toy model of a single bosonic degree of freedom with derivative self-interactions. 

We consider the real-time evolution of a condensate and derive evaporation laws for two possible 
decay mechanisms in the Schwinger-Keldysh formalism. We show that semiclassical results can be 
reproduced if the decay is due to an effective two-body process, while the existence of a three-body 
channel would imply very short lifetimes for the condensate. In either case, we uncover the existence 
of scaling solutions in which the condensate is at a critical point throughout the collapse. In the 
case of a two-body decay we moreover discover solutions that exhibit the kind of instability that 
was recently conjectured to be responsible for fast scrambling. 


I. INTRODUCTION 

Black hole (BH) solutions in general relativity (GR) do 
not possess global timelike Killing vector fields. As a con¬ 
sequence, black holes of mass M emit radiation [I], semi- 
classically found to obey a thermal^ distribution with a 
temperature Th (x 1/M. In a semiclassical analysis, an 
evaporating black hole can turn pure states into mixed 
states [2] - the evaporation process appears to violate the 
unitarity principle of quantum mechanics. 

More recent advances, most prominently in the holo¬ 
graphic [3] understanding of black holes have led to the 
widespread belief that black hole evaporation is after all 
a unitary process. In the above terms, this necessarily 
implies a breakdown of the semiclassical description at 
some point during the evaporation process, most likely 
around the halfway point of evaporation when the black 
hole is still large [4]. In particular, this calls for an exten¬ 
sion of the semiclassical approximation towards a more 
complete picture. 

In this paper we develop a microscopic mechanism re¬ 
sponsible for particle loss of black holes within the con¬ 
densate picture that has been proposed in [5, 6]. There, 
black holes are described as multiparticle quantum states 
of gravitons at a point of collective strong coupling, or 
quantum criticality. The key assumption is that all rel¬ 
evant physics originates in collective effects of gravitons 
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^ Of course, due to the decrease of the black hole mass, the ob¬ 
served spectrum will not be thermal for a finite mass black hole. 
Exact thermality is a consequence of the semiclassical limit cor¬ 
responding to black hole mass M —> oo, Planck mass Mp —> oo, 
and Schwarzschild radius Rs = M/Mp fixed. 


whose wavelength is given by the characteristic length 
scale t of the gravitational background. This in principle 
allows a qua ntum description of objects of size t Ip, 
with ip = 'JhG the Planck length, within the low en¬ 
ergy approximation to the effective action of gravity, i.e. 
the Einstein-Hilbert action. Within this description, two 
effects can lead to significant deviations from semiclas¬ 
sical predictions. First, the criticality of the condensate 
may ultimately lead to the appearance of almost gap¬ 
less modes, corresponding to collective excitations of the 
graviton condensate. Their presence can lead to strong 
quantumness even for large black holes [7]. Second, cor¬ 
rections suppressed by the effective number of gravitons 
may become important on sufficiently long timescales. 
This, in turn, serves to restore unitarity in the evapora¬ 
tion process [6]. 

In the case of black holes, wavelengths of the order of 
the Schwarzschild radius Rg are expected to dominate. 
Under this assumption, black holes indeed lie at a point of 
collective strong coupling. This may be the microscopic 
origin of black hole entropy [6] and the underlying reason 
for black hole quantumness [7] and scrambling [8]. 

The strength of graviton-graviton interactions is given 
by perturbative general relativity as 


where £ is the characteristic wavelength of the interact¬ 
ing gravitons. For large systems, the particles are ex¬ 
tremely weakly coupled. However, due to their bosonic 
nature, they can occupy states with a large number of 
particles. In a Hartree picture, each graviton sees a col¬ 
lective binding potential produced by the other gravi¬ 
tons. For sufficiently large numbers, this allows for the 
formation of self-sustained bound states. The number of 
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constituents N of such a bound state can be estimated 
using the virial theorem [5], where self-sustainability is 
achieved for aN = 1. The average wavelength of a gravi¬ 
ton in the system is related to N via Eq. (1) 

£ = ip^/N, ( 2 ) 

while the total mass of the bound state is approximately 
given by the sum of the energies of the individual con¬ 
stituents M = Nj = y/N As a consequence, N be¬ 
comes the universal characteristic of the condensate. 

Black hole formation from gravitational collapse may 
be understood as bringing the condensate to the critical 
point [6]. In fact, N is set by the energy stored in the 
gravitational field, roughly constant throughout the col¬ 
lapse. On the other hand a increases until black hole 
formation sets in for aN ~ 1. The formation of BH 
in ultrahigh energy scattering at center of mass energies 
-y/s ^ Mp also appears accessible within this approach. 
Graviton numbers of order l/a ~ s/Mp are seen to dom¬ 
inate scattering amplitudes, either as intermediate states 
in 2 —> 2 processes [9], or as final states in 2 —)■ TV scat¬ 
tering [10]. 

Condensate decay, on the other hand, is due to the 
interactions between the constituents. In essence, inco¬ 
herent scattering between condensed gravitons leads to 
ejection from the condensate and a change of particle 
number N ^ N'. Concurrent to the emission, the black 
hole collapses and readjusts its size: —)■ i?[, = IpVW. 

Consequently, the graviton-graviton coupling is given by 
a = ipjR's = l/N'. After one emission cycle, the con¬ 
densate is still at a point of collective strong coupling 
with lower particle number and stronger coupling. 

In this paper we address the mechanism for the evapo¬ 
ration of the condensate in a microscopic picture (see e.g. 
[11, 12] for other considerations towards Hawking evapo¬ 
ration in this framework). We consider a toy model with 
an interaction structure that bears important similarities 
with GR. We demonstrate that the evaporation of the 
condensate is to leading order due to scattering of two 
condensed particles. We further show that if one of the 
two participants can rescatter into the condensate, the 
obtained rate greatly exceeds the rate that is expected 
from semiclassical considerations. If, on the other hand, 
this is disallowed, as for example for a homogeneous con¬ 
densate, we obtain an evaporation rate that scales as 

or, using M ~ MpV N, in terms of the bound state mass 

M = -j^+OiN-^). (4) 

In the large N limit, this corresponds to the mass loss 
that is obtained semiclassically for a black hole of tem¬ 
perature Th = yR/n picture is applicable to 

gravity, unitarity is obviously never an issue; the entire 


black hole can radiate away without any loss of informa¬ 
tion. The difference to semiclassical results is encoded in 
1/A^ corrections. 

Moreover, we discover that in the case of a two-body 
decay, the critical solutions exhibit an instability with 
Lyapunov coefficient A ~ 1/lp^/N. This implies that 
along these solutions, entanglement is generated on a 
timescale Rs log i?s. It has been conjectured [6] that this 
is the microscopic origin for the fast scrambling property 
of black holes [13]. 

Our paper is organized as follows. We discuss our as¬ 
sumptions in section II. In sections III and IV we de¬ 
scribe the basic physical mechanism responsible for the 
evaporation of the condensate. We derive the equations 
governing the collapse of the condensate in section V A 
and take evaporation into account in section VB. The 
final section then contains discussions on the solutions 
found in case of three- (VIA) and two-body (VIB) de¬ 
cay. From here on, the speed of light c and Planck’s 
constant h are chosen to be unity. 

II. FROM GRAVITY TO PROTOTYPE 

In the condensate picture for black holes, the dynam¬ 
ics of collapse and Hawking evaporation are due to two 
intertwined effects. The coherent excitation of a tachy- 
onic breathing mode of the condensate leads to collapse 
of the black hole. This is a process involving only the 
gravitons of the condensate. At the same time, incoher¬ 
ent scattering allows for the production of gravitons that 
can escape the black hole. In principle, the former may 
be accounted for through mean field evolution, while the 
latter is due to the interaction of the mean field with 
quantum fluctuations. 

In a gauge where the linear^ graviton obeys the rela¬ 
tions 

hoi = 0, /i]) = 0, = 0, (5) 

the corresponding time evolution is generated by a 
Hamiltonian that in Fourier space takes the (suggestive) 
form 

H= fd^kJ2 |k|<Aak,A±^M2-"x 

A^l,2 n 

f d^ki... cPK , (g) kA . (6) 

Here A = 1, 2 corresponds to the two transverse polariza¬ 
tions of the graviton and the at ,a are creation and 

annihilation operators of gravitons with polarization A. 
The functions comprise all possible degree n mono¬ 
mials of said creation and annihilation operators, thereby 


^ As usual, we linearize around Minkowski. This also gives us a 
preferred time slicing. 
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generating the infinite series of vertices present in an in¬ 
teracting massless spin-2 theory. The interaction term 
will generically contain also the longitudinal and tempo¬ 
ral polarizations of the graviton, depending on the choice 
of nonlinear gauge. 

In our picture, a black hole roughly corresponds to a 
set of quantum states |BH) in the interacting theory with 
a large occupation number of gravitons in a single mode 
Q-BH = / d^kauak- Note that the dfc comprise annihila¬ 
tion operators of all possible polarizations, in order for 
the black hole state to be part of the physical spectrum. 

In principle, (6) contains all the relevant information 
for the analysis of the non-equilibrium behavior of the 
states |BH). However, the corresponding vertices, stem¬ 
ming from a Poincare invariant field theory, will not con¬ 
serve particle number. Moreover, we are dealing with the 
presence of the infinite series of vertices and the ambigu¬ 
ity due to gauge redundancy. An explicit treatment of 
this theory is extremely difficult. Therefore, we are in 
dire need of simplifying assumptions. 

For now, these will be: 

(i) Reduce the number of polarizations to a single 
mode, thereby also removing the gauge ambiguity. 

(ii) Focus on particle number conserving processes. 

(iii) Assume that the relevant dynamics of the conden¬ 
sate is captured already by the lowest order inter¬ 
action. Due to (ii) we consider only the quartic 
vertex. 

We will briefly comment on their viability in the Ap¬ 
pendix. 

The latter two assumptions will prevent us from learn¬ 
ing anything about the actual spectrum of emitted parti¬ 
cles. Nevertheless, we will see that they already allow for 
very interesting conclusions on the condensate dynamics. 
Take note here that due to (ii) and (iii), all momenta 
involved in the leading order collision processes are ex¬ 
pected to be of the same order. We will therefore also 
neglect the momentum dependent prefactor of the quar¬ 
tic interaction term. 

With these assumptions, we arrive at the Hamiltonian 



- J d^ki .. (X! 

Under our assumptions, the difference to previous pro¬ 
totype models for graviton condensates [6-8] reduces to 
the relativistic dispersion relation. We will see, however, 
that it is precisely this feature that is responsible for in¬ 
teresting properties. Note also that our Hamiltonian is 
based on derivative interactions. The difference in the 
interaction as compared to [14] is due to the inherently 
nonrelativistic nature of the model considered there. 

Before we start analyzing the dynamics of (7), we utter 
a word of caution. Of course, the simplified Hamiltonian 


is void of quite a few important features of gravitation. 
Besides the simplifications in terms of number conser¬ 
vation, we have eliminated the longitudinal and tempo¬ 
ral modes from the dynamics. In GR, it is precisely 
these modes that are responsible for the gravitational 
potential. Most simply, their presence may be modeled 
through the inclusion of a trapping potential. This will 
be left for future work. Let us also note that the reduc¬ 
tion to particle number conserving processes simplifies 
the structure of the interacting vacuum. This will al¬ 
low us for example to decompose correlators in a coher¬ 
ent state basis without having to worry about subtleties. 
Implications of the nontrivial nature of the interacting 
vacuum on the bound state description of black holes 
have been studied in [15]. 


III. SCHWINGER-KELDYSH EORMALISM 


The following section provides an elementary review of 
the Schwinger-Keldysh [16] formalism that will allow for 
a proper treatment of the real time dynamics of collapse 
and evaporation. It may be skipped by the experienced 
reader. Comprehensive introductions may be found e.g. 
in [17]. Our analysis presents a generalization of previ¬ 
ous results (found e.g. in [18-20] and many subsequent 
works) to systems with a relativistic dispersion relation. 

We are interested in the time evolution of an unstable 
condensate which is initially^ described by a (normal¬ 
ized) density matrix p{ti). The expectation value of any 
observable O is given by 


(0)(t)=Tr 


u{u,t)du{t,ti)p 


= Tr 


u{u,tf)u{tf,t)du{t,ti)p , 


( 8 ) 

(9) 


where U{ti,t 2 ) = exp — ^ 2 )) is the time evolu¬ 

tion operator. The second equation has been obtained 
through an insertion of U{ti,tf)U{tf,ti) for tj in the 
asymptotic future. It has served to extend the integra¬ 
tion path from ti ^ t ^ ti to ti ^ tf ^ ti, the so-called 
Keldysh contour, which we will denote by C. 

With use of the Keldysh contour, expectation val¬ 
ues can be obtained from a generating functional 
through the introduction of corresponding sources into 
the Hamiltonian 


Z[J] 

{0){t) 


H±J{t)0, 

(10) 

Tr[[/j(C)/5] , 

(11) 

1 6Z[J] 

2 

(12) 


^ A precise study of the dynamics of black hole formation in the 
condensate picture is an interesting issue on its own. Here we 
simply assume the presence of an initial condensate. 
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Here and H~ are the Hamilton operators along the 
forward and backward contour, respectively. 

As usual, (11) may be turned into a path integral by 
introducing at each timestep an appropriate partition 
of unity. In this case, we use coherent states that are 
eigenstates of the annihilation operators appearing in (7), 

Ok IV') = V'k IV’)- One obtains 

Z = j , (13) 

where the “action” S is given by 

S'[V’,V’^] ^ J y c?^fci|#kj'9tV’ki - IkilV’kjV’ki 
+ J Sk2 ■.. d^fc4V’Ii V’l2V'k3V’k4<5^^^ (X! ^0 I • 

(14) 

The information on the initial state is encoded in the cor¬ 
relation of the field on the forward and backward branch. 

The time integral in (14) may be brought into conven¬ 
tional form by introducing forward and backward fields 
V’^ that live on the forward (backward) branch of the 
Keldysh contour. Performing a so called Keldysh rota¬ 
tion by introducing the “classical” and “quantum” fields 

V>'=^ = y= (V^+-hV'") , V’'^ = (V’"^ - V’") , (15) 

one obtains 

J dt d^kid^k2|iv)k^/c(ki, k2)v)k2 
+ Mp-2 J ki) 

X (V’k?V’kJV’k3V’k4 + V’klV’kIV’k3V’£?) I + h-c- • (16) 

Here we have introduced V’k = (V’k > V’k)^ ll^® kinetic 
matrix 1C is defined as 

/C(ki,k2) = 

0 d(ki-k2)(fat - |ki|)\ , . 

d(ki-k2)(z9t-|ki|) S;^(ki,k2) )■ 

We have introduced the Keldysh self-energy S/f, whose 
precise value depends on the interactions and is of no 
particular interest to us. Its presence, however, is im¬ 
portant, since it contains the information on the correla¬ 
tors of the forward and backward fields and thereby on 
the initial density matrix. Note that at this point, one 
may equivalently seek a formulation for the action (16) 
in terms of real fields by translating 4' and 4*1 into a real 


scalar field and its canonical momentum. However, our 
focus on particle number conserving processes is more 
straightforwardly implemented in the current language. 

The classical mean-field dynamics of the condensate 
arise from (16) as the solution to the saddle-point equa¬ 
tions that has V’'^ = 0- In this case, variation of (16) 
with respect to V’'^^ yields the Gross-Pitaevskii equation 
for the classical field. This is the equation of motion that 
describes the mean-field dynamics of the condensate: 

= |k|V’£' 

- M-2 j d^k.d^k^d^ksS (k + ki - k2 - ka) 

( 18 ) 

For a condensate with attractive self-interactions, the 
normalized solutions to this equation will correspond to 
a collapsing condensate of N particles once N surpasses 
a critical value [21]"^. 

IV. EFFECTIVE ACTION 

The evaporation of the condensate is due to scatter¬ 
ing of condensed particles into the quasi-particle cloud. 
These effects may be taken into account most readily by 
integrating out the quasi-particles. To this regard, we 
separate the fields into a condensate part and fluctua¬ 
tions. 

= (j> + 6ip, + S(p^ . (19) 

Inserting (19) into (16) gives rise to a plethora of terms, 
of which only few have a relevant effect. The reason for 
this lies in the fact that we are dealing with a condensate 
with V 1 at T = 0. The quasiparticle occupation is 
much lower than that of the condensate mode. We may 
therefore focus on terms that are at most quadratic in the 
fluctuations. All corrections that arise from higher order 
terms are proportional to the density of fluctuations and 
are thus at least 1/V-suppressed. 

Primarily, we are interested in the loss of condensed 
particles due to incoherent scattering. Two processes 
contribute: 

(i) Two-body decay: This corresponds to scattering of 
two condensed particles and subsequent emission 
of two quasiparticles. In the action the relevant 
operator is quadratic in the fluctuations. 

(ii) Three-body decay: Again, two condensed particles 
scatter. Here, however, only one of the outgoing 
particles leaves the condensate. The corresponding 
operator is linear in the fluctuations. 


^ Even for smaller N, the condensate is only metastable and can 
always decay through macroscopic tunneling [22]. 
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Note that the three-body process is only non-vanishing 
for an inhomogeneous condensate; otherwise it is incom¬ 


patible with the conservation of momentum. 

Following our above remarks, we construct the follow¬ 
ing expression for the quasiparticle action: 


5qp = 


dtz(5^k^/Cqp(ki,k2)d0k2 + 


-2 


d^k; 


id^k 




X 


(26(j)' 


cli 






■ 6(j)] 


ki '/’L ^ks ^ki 




^L4.</>k3 



-I- h.c., 


( 20 ) 


We have here introduced the quasiparticle vector 5(j) = 
(^S(I)^\S4>'^, as well as the quadratic quasi¬ 
particle operator /Cqp, which reads 


fJC +A B \ 

\ K."! + a) ^ 


( 21 ) 


with fC defined in Eq. (17) and (suppressing momentum 
labels and integrals) 


A = 2Mp-2 






B = Mp-2 






-f h.c., 


( 22 ) 

(23) 


We have retained only terms that are at maximum lin¬ 
ear in the quantum field ())*. A loss term in the Gross- 
Pitaevskii equation can only originate from a term in the 
effective action that is linear in all other effects that 
may alter the condensate dynamics are 1/A^ suppressed. 
Note that the quadratic part of the fluctuation action 
allows one to read off the quasiparticle spectrum; diago- 
nalization leads to the celebrated Bogoliubov modes [23]. 

In order to obtain a modified Gross-Pitaevskii equa¬ 
tion that incorporates the effects of the quasiparticles, 
we integrate out the latter. Loss terms will be generated 
by the diagrams shown in Fig. 1, which obtain imaginary 
parts due to on-shell fluctuations. Fig. la describes the 
three-body decay process [19], while Fig. lb corresponds 
to the two-body decay. In the effective action, we gener¬ 
ate the contribution 


^^^'kid3k2r(ki,k2) - KA.) ’ 

(24) 

where we have introduced 


r(ki,k2) = r( 3 )(ki,k 2 ) + r(2)(ki,k2) (25) 


and defined 

r(3)(ki,k2) = 

J d^M^ksd^kiS {\ki\ + \k2\ - [k] - |ki +k2-k|) 

^ '^ki-|-k3-k')’k3 0k^0k2+k4-k , (26) 

r( 2 )(ki,k 2 ) = 

J d^kd^kgJdkil + |k2l - [k] - |ki + k2 - k|) 

^ ‘t^ki+k^-k^^ks ■ (27) 

Here the former contribution is due to diagram la, giving 
rise to a loss term in the effective action that is propor¬ 
tional to six powers in the fields, hence scaling with the 
particle number as N'^. On the other hand, the contri¬ 
bution from Fig. lb is proportional to only four pow¬ 
ers in the fields, scaling as . Note that expression 
(27) is obtained from Fig. lb by evaluating the integral 
/ dB^5{k^ — jk|), following from the on shell contribution 
from the propagator of k. Let us here also mention that 
one-loop diagrams with more external legs correspond 
to processes in which one or more of the emitted quasi¬ 
particles scatter on the condensate. These give rise to 
contributions which are proportional to the quasiparticle 
density; they are thus 1/Wsuppressed and will be ne¬ 
glected. Such additional suppression is also present for 
all higher loop processes. 

In order to ensure that we do not erroneously count 
processes where all particles rescatter into the conden¬ 
sate, we should impose a constraint on the fluctuation 
momentum k and only integrate over momenta outside 
the support of the condensate. This is vital in particu¬ 
lar when considering a homogeneous condensate, where 
this ensures that the three-body decay rate vanishes. On 
the other hand, for an inhomogeneous condensate, this 
will only lead to a numerical renormalization of the decay 
rate, unless the window of allowed momenta depends ex¬ 
plicitly on N. In absence of a compelling argument why 
this should be the case, there will always be a critical N 
after which the three-body decay dominates. 

Once the above loss terms are included, the Gross- 
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(a) Three-body decay 




(b) Two-body decay 


FIG. 1: Lowest order diagrams responsible for the imaginary part of the self-energy. The external lines correspond 
to condensed particles while internal lines correspond to quasiparticles; solid lines identify classical and dashed lines 
quantum fields. The imaginary parts are induced when in (a) k and in (b) k and fci + ^2 — fc go on-shell. 


Pitaevskii equation reads 
= |k|(/)k 

+ J d^kid^k2fi^k3 (k + ki - k2 - ka) 0k^(/)k20k3 

-i J (i^k2r(k, k2)(/)k2 • (28) 

From this, we can immediately read off the change in the 
number of condensed particles 

^ = -J d'kid3k2F(ki,k2)<^0k2 (29) 


variational solution to the simpler problem without the 
dissipative term. This will provide us with an equation 
for the condensate size. Second, we supplement this with 
Eq. (29) in order to take into account the loss of particles. 
There we will focus on two regimes. First, we consider 
the case in which the three-body decay (26) is allowed. In 
the second case, the dominant contribution comes from 
the two-body decay (27). 

V. VARIATIONAL APPROACH 
A. Collapse 


Having integrated out the quasiparticles, the final step 
towards an effective action that describes the collapse 
and evaporation of our condensate is now to integrate 
out the quantum field (/)*. We can do this via integrating 
in an auxiliary “noise” field ?7k(0- 

Z = j V(j) V(j)^ Vr] , (30) 

Sn[(l),(l)^] = J - ryk) +h.c. 

+ (31) 

where £ is the operator corresponding to the Gross- 
Pitaevskii equation (28). Integrating out (j/^ now con¬ 
strains the classical field (j)‘^ to obey a Langevin equation 
with Gaussian noise rj 

£wi(p,<l)^) =Vk- (32) 

The dynamics described by (32) may in principle be 
obtained numerically. Instead, we will here take a sim¬ 
pler route and seek a variational solution to the aver¬ 
aged Gross-Pitaevskii equation. Taking note that the 
dissipative equation (32) does not directly follow from a 
variational principle, we proceed by dividing the varia¬ 
tional approach into two steps: First, we shall look for a 


In the following we will focus on the average dynamics 
of the condensate. By the Gaussian nature of the noise 
ry, its first moment vanishes, and will hence not influence 
the dynamics. Without the dissipative term, Eq. (32) 
can be obtained as the stationary point of the following 
Lagrangian. 


/Ik = ^ - |k||0kP 

+ M-^J d3kid3k2d3k3<5(3) <(-k2<)>k3) 

+ h.c., (33) 


We extremize Eq. (33) with respect to a set of spheri¬ 
cally symmetric trial functions. Guided by simplicity, we 
choose a Gaussian ansatz in real space^ 


<p{r,t) = A{t) 



e 


2L(3Y( -ri- 
2 UJ \R(t)' 



(34) 


A{t) is the complex amplitude, R{t) the real width of the 
condensate and r the radial coordinate. The function h(t) 


^ Note that for a nonrelativistic harmonically trapped condensate, 
the ground state wavefunction can indeed be well approximated 
by a Gaussian even in the presence of interactions [24, 25]. 




7 


can later be identified with the velocity of the collapse. 
The normalization of (f) is chosen such that = N 

and the numerical factors simplify the calculation. 

Fourier transforming Eq. (34) and inserting the result 
into Eq. (33), we obtain the averaged Lagrangian density 
L = / d^kC 


L = - AA*) - ^\A\'^R%- 


2R 


\/l + IPR^ + 


27|A|^ 


128V2M^R3 ■ 


(35) 


We can now understand the collapse dynamics as a 
variational problem of the time dependent parameters 
s = {A, A*, R, b}, which obey the equations of motion 


dt \ dsi) 



(36) 


The equations of motion for the amplitude simplifies 
to particle number conservation 

^|Ap = 0, (37) 

or in other words N = 0. 

Next, we can relate b to the collapse velocity R by 
varying for b and using Eq. (37) 


b = 


R 


(38) 


This equation completely determines the function b in 
terms of R and R. 

The expression for the condensate width is obtained 
by variation with respect to R. After substitution of 
Eq. (37) and (38), one obtains 






V2 V4 


N 


M2 


(39) 


Note the Lagrangian after integrating out b 


2R 72 




(40) 


We recognize it as the Lagrangian of relativistic point 
particle with a mass that depends on R. 


1. Slow collapse 

It is instructive to first consider the case of small col¬ 
lapse velocities, as this corresponds to a nonrelativistic 
limit and allows us to qualitatively compare our expres¬ 
sions with existing results from the literature. 



In the limit of small velocities and small acceleration, 
Eq. (39) reads 

- 1 72 I fsV N 

R^ 2R 72 Uy ■ 

The first term in Eq. (41) corresponds to the outward 
force due to the kinetic energy of the bosons; its scaling 
with inverse R is dictated by Heisenberg’s uncertainty 
principle. The third term is due to the attractive inter¬ 
actions and can, for sufficiently large N, overcome the 
repulsive force. In that case, the condensate collapses. 
In comparison with known results in the literature (e.g. 
[19, 24]), the second term corresponds to corrections due 
to the relativistic dispersion relation. 

In the small velocity limit, the Lagrangian (40) may 
be canonically normalized. The corresponding equation 
of motion will then give us a simple picture of the time 
evolution of the width as the motion of a particle mR = 
— -j^V{R, N) in a one dimensional potential. To see this, 
let us take the small velocity limit also in the Lagrangian 
(40): 


^_3^^72 3N 27 

12872 M2i?3 ■ 


(42) 


We may canonically normalize the kinetic term through 
the redefinition R = R^. From the canonical Lagrangian 


L = 3IVI?2 


3N 1 /3y A^2 

^ 7^ Uy Mp6 ’ 


(43) 


we can conclude that the motion corresponds to that of 
a particle with mass m = 6N in the effective potential 


V{R) = 


3N 


1 /3y 

71 Uy Mp® ■ 


(44) 


We plot the effective potential in Fig. 2. The potential 
V{R) possesses a maximum, located at 


Rl =R+ r^lpy/N. 


(45) 

















It turns out that for R ~ i?_|_, the criticality condition is 
fulfilled: 

/2 

aN ~ -^N - 1. (46) 

Finally note that the runaway behavior for large R is 
due to the fact that we have not included an external 
trapping potential for the condensate. 


B. Evaporation 

The decay rate of the condensate according to Eq. (29) 
can be evaluated for a spherically symmetric collapse us¬ 
ing the ansatz Eq. (34) which yields 

Here the C(i) are dimensionless constants, whose value 
depends on the precise form of the condensate wavefunc- 
tion. The factors of Mp can be read off straightforwardly 
from the diagrams Fig. 1, while the dependence on R 
follows from dimensional grounds. 

As we have discussed, in our setup both C( 3 ) and C( 2 ) 
are nonzero. In that case, for sufficiently large N the 
evaporation will be dominated by the three-body decay 
with N ~ —N^. However, we admit the possibility that 
the precise form of the interaction or the form of the 
wavefunction can effectively close the three-body decay 
channel. In this case, the evaporation will occur with 
N ~ —N"^. Both cases exhibit very interesting proper¬ 
ties, which we shall explore in the next sections. 

Equations (39) and (47) dictate the evolution of the 
condensate described by the Hamiltonian (7) in the Gaus¬ 
sian approximation, as long as higher order correlators of 
the fluctuations can be neglected. 


VI. SOLUTIONS 
A. Three body decay 


Two properties of these solutions are immediately evi¬ 
dent. Eor one, they exhibit self-similarity: The critical¬ 
ity condition R ~ ipy/N is fulfilled throughout the col¬ 
lapse. Second, the total evaporation time is proportional 
to tcoii ~ ~ £p\/Nq. Irrespective of the collapse ve¬ 

locity u, this immediately rules out a three-body decay 
of the form (26) as the mechanism of black hole evapora¬ 
tion. In gravity semiclassical arguments indicate a black 
hole lifetime teH ~ ~ which is a factor 

Nq longer than the lifetime of the condensate considered 
here. 

Let us nevertheless comment on some of the interesting 
features of the solutions (48). We find three solutions in 
the allowed range 0 < u < I; a c-independent solution 
u = I as well as the two c-dependent ones®. Let us point 
to the curiosity that the two latter solutions are only real 
for values of C( 3 ) that are smaller than a critical value 
Ccrit- At C( 3 ) = Ccrit both solutions disappear in a saddle- 
node bifurcation. We have illustrated this behavior in 
Eig. 3a. 

Let us now investigate the stability properties of the 
critical solutions (48). As we will see, not all three solu¬ 
tions are stable. However, as long as the two c-dependent 
solutions exist, at least one of them presents an attractor. 

Since we are only interested in deviations from critical 
behavior, and not in the stability of a specific solution, 
it is sufficient to decompose R as 

R{t) = ( 1 ^) vW)tp + mt) , ( 50 ) 

and to then combine Eqs. (39) and (47) and linearize in 
5R. Eor V = I, we obtain 


5R = 


£IN 


2 

C(3) 


6R -= 

£pVN 


2 

C(3) 


SR, (51) 


while in the other cases, we have 

12r;2 _ 2 


dR = 


£IN 


^Ysr- 

C(3)/ 


5R. 

£p\/N Vc(3) / 


(52) 


One may solve the collapse and evaporation equations 
(39) and (47) numerically for generic initial conditions, 
drawing a complete picture of the behavior of the con¬ 
densate in the variational approach. However, it turns 
out that in case of the three-body decay, the equations 
possess a simple set of analytic collapse solutions: 


The eigenvalues of the stability matrix read 



-3 ± V 29 

23/4c(3)1/4£pV^ 


in the first case and 


R{t) = Ri — vt 


Nit) 


I 2v (Ri — vt)'^ 

\] <3(3) £p 


(48a) 

(48b) 


A 


(2.3) 

± 


2U2.3 

C(3) 


-3^2,3 ± ^2 - 3u|3 

£pVN 


(53) 


(54) 


The parameter v is fixed via the algebraic relation 



® The corresponding values are given by n = cos (/(c)) di 
sin (/(c)), where /(c) = |arctan (^/(|)^ ^ 



9 


V KeX+£p^/l^ 




(b) Largest eigenvalue as a function of c. The solution u = 1 (blue, 

(a) Velocity parameter u as a function of c. While the solution u = 1 long-dashed) always possesses an unstable direction; the solution 
(blue, long-dashed) is c-independent, the solutions with v = U 2,3 v = V 2 (black, solid) is absolutely stable, while v = V 3 (red, 

(black, solid and red, dashed) are only real for c < Ccr- For c = Ccr, short-dashed) is unstable for small c. For c < 0.05, are 
the solutions disappear in a saddle-node bifurcation. degenerate. This is the origin of the kink in the black curve. 

FIG. 3: Solution properties for the three-body decay. 



FIG. 4: Numerical results for R(N) for initial conditions 
close to the critical point i?in ~ \/Nin and c = O.OI. 


in the latter. The expressions on the corresponding so¬ 
lutions for V are rather lengthy. Important here is only 
that real parts of the eigenvalues are always negative for 
V 2 , while for V 3 the larger one is positive for c below 
some threshold value. Henceforth, at least one of the 
two solutions is absolutely stable. The solution = 1, 
on the other hand, possesses an unstable direction; un¬ 
der small perturbations, it flows towards the solution 
n « 1 — 2.8 c(3 ). The scaling behavior therefore remains 
unaltered. We display the stability behavior in Fig. 3b 
by plotting the real part of the larger eigenvalue for all 
three solutions. 

To complete the picture, we present numerical results 
for initial conditions close to criticality in Fig. 4. As 
expected, we observe that all solutions flow towards the 
critical solution with q{N) ~ '/N. 

Their peculiar properties in terms of stability and crit¬ 
icality make condensates that behave according to (48) 
interesting objects on their own. Ultimately, however. 


we are interested in a condensate that could mimic the 
behavior of black holes. Their very short collapse times 
rule out condensates that deay via a three-body process. 
We therefore now focus on condensates whose dominant 
decay mechanism is given by the two-body process Fig. 
lb. 


B. Two body decay 

In the case of C( 3 ) = 0, the simple solutions (48) are ab¬ 
sent. Nevertheless, we find solutions that approximately 
display critical behavior. 

In order to see this, we write R as 

R = Rc{N) . (55) 

Plugging this into Eqs. (39,47) with C( 3 ) = 0 and de¬ 
manding mutual consistency of the equations yields a 
differential equation for Rc{N). In absence of an exact 
solution to this differential equation, we write the rela¬ 
tion between N and i? as a power series: 

R,{N) = ip^/^ ^ • ( 56 ) 

Expanding in powers of 1/N, we hnd up to second order 

RciN) =(^l- (57) 

with qq = (1)^- The criticality condition is hence 

fulfilled up to corrections of order 1/A^^. 
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To lowest order in 1/A^, we find the solution for the 
particle number 

2 

jV(f)^iVo (l-| 3 I +0{\/N). (58) 

Not surprisingly, if we were to trust this solution up to 
the point of complete evaporation, we would obtain an 
evaporation time t ~ ~ ji^p. 

Again, we test the stability of the criticality condition 
by considering 


R{t) = Rc{N{t)) + 6R{t) , (59) 

linearizing in SR and including contributions up to order^ 
1/N'^. We obtain the two eigenvalues 





2al^^£pN^/'^ ' 


(60) 


For large N, A+ is positive. The critical solutions are 
thus unstable. The leading Lyapunov coefficient is found 
to be Xl ~ 1 /ipVN. 

This result opens up a curious side track. As was 
shown in [8], the existence of an instability in the Gross- 
Pitaevskii equation leads to generation of one-particle en¬ 
tanglement on a time scale t ^ Xj^^logN. Using the 
above result, we obtain t ^ i?logi?, implying that quan¬ 
tum correlations will become important on timescales of 
the order of the scrambling time of black holes. This 
could provide a first hint towards a fast scrambling be¬ 
havior of these kind of condensates. 

In the case of black holes, the solution should follow 
the critical behavior at least up to the point in which 1 /N 
effects can become dominant, a time roughly of the order 
of Page’s time [4]. In the language of gravitons, this is 
due to the maximal packing property [5] of BH, which 
reflects the fact that it is impossible to localize energy in 
a volume smaller than its own Schwarzschild radius. This 
kind of behavior cannot be concluded from our analysis. 
This may be due to several reasons. For one, it may be 
that maximal packing is intimately related with entropy. 
In the condensate picture, exponential degeneracy is only 
expected to be present at the critical point. This could 
imply that trajectories that follow critical behavior are 
quantum mechanically favored. If the large entropy is a 
unique feature of GR it can never be seen in simpler toy 
models. This may be tested through an analysis of “clas- 
sicalizing” [26] theories, which are expected to show simi¬ 
lar behavior. Most probably, this would require applying 
our analysis to a theory with more powers of derivatives. 

On the other hand, it is conceivable that the apparent 
absence of maximal packing is not due to the toy model. 



FIG. 5: Numerical results for R{N) for initial 
conditions close to the critical point Rq ^ The 

deviation from critical behavior becomes relevant at 
times t ~ i?o log Ro- 


but instead due to the approximation in which we ne¬ 
glect higher order correlators in the Bogoliubov hierar¬ 
chy. In our toy model, the instability implies a deviation 
from critical behavior that becomes 0 { 1 ) on the same 
timescale as the deviation from classicality. On the other 
hand, this in turn implies that corrections to Eqs. (39,47) 
can become relevant. A conclusive answer requires fur¬ 
ther analysis. 

For completeness, we show in Fig. 5 the behavior of 
the condensate for generic initial conditions in the vicin¬ 
ity of the critical solution. Indeed, we observe that 
after a time t ~ R\ogR, the solutions diverge. For 
R{tin) > Rc{N{tin)) the condensate eventually ceases to 
collapse and turns around, while particle emission con¬ 
tinues for a while. At some point, the condensate is so 
dilute that emission is effectively shut off and the par¬ 
ticle number remains constant. On the other hand, for 
R{tin) < Rc{N{tin)), we observe a rapidly collapsing con¬ 
densate after the instability time. The emission remains 
close to the solution (58) until the condensate size has 
reduced to a fraction of its initial value. Then, almost all 
condensed particles are ejected within a very short time. 

However, we stress here once again that due to the 
presence of the instability, we expect higher order corre¬ 
lations to become large around the instability time. The 
description in terms of an effective Gross-Pitaevskii equa¬ 
tion is then likely to break down®. 


VII. CONCLUSIONS 

In this work, we have developed a toy model for Hawk¬ 
ing evaporation in the context of the Bose condensate 


^ Note that if we expanded to higher orders, source terms would 
appear that are due to (56) deviating from the exact solution to 
the equations of motion. 


Taking into account all corrections would lead to the well known 
BBGKY-hierarchy [27]. In practice, the hierarchy has to be trun¬ 
cated at some order; in the case of rapidly growing correlations, 
this is difficult to do consistently. However, this is certainly a 
possible road that deserves further attention. 
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picture for black holes. To this end, we have constructed 
a Hamiltonian that captures essential ingredients to the 
underlying physics of black hole evaporation, while at the 
same time being stripped down from some of the com¬ 
plications that arise in Einstein gravity. In particular, 
we have focused on a single degree of freedom and have 
turned off processes that violate particle number conser¬ 
vation. We have introduced the Schwinger-Keldysh for¬ 
malism for nonequilibrium dynamics of Bose condensates 
and derived the Keldysh action that describes a collaps¬ 
ing and evaporating condensate. 

We have then chosen a variational approach to solve 
the ensuing equations of motion, using the number of 
condensed particles N and the size of the condensate R 
as variational parameters. The resultant action took the 
form of that of a relativistic particle in an effective N- 
dependent potential. In comparison with existing results 
in the literature, we have identified the corrections due 
to the relativistic dispersion relation. 

Once the decay of the condensate was taken into ac¬ 
count, we have found two possible behaviors. If the con¬ 
densate is sufficiently inhomogeneous, evaporation is pre¬ 
dominantly due to interactions in which one of the par¬ 
ticipants rescatters into the condensate. In this case, the 
resulting decay rate is proportional to Vn jIp ^ leading to 
a lifetime of the bound state that is much shorter than 
what is expected for a black hole. Within the regime 
of validity of our assumptions, the three-body decay is 
required to be forbidden if the condensate should be a 
viable black hole candidate. 

On the other hand, we have considered situations in 
which the three-body process is forbidden, for example 
by momentum conservation for a homogenous conden¬ 
sate, and decay can only occur through a two-body pro¬ 
cess in which both particles are ejected. In this case, 
we have discovered many of the properties that make a 
condensate picture for black holes appealing. Collapse 
and evaporation can happen nearly self-similarly, with 
width and particle number related via R ~ lp\fN up 
to subleading I/iV corrections. The bound state then 
has a lifetime that may be as long as much like 

semiclassical black holes. At the same time, the solution 
exhibits an instability of the form conjectured in [8] to be 
responsible for scrambling. The leading Lyapunov coeffi¬ 
cient was found to be l/tp\/N, yielding a quantum break 
time of order R log R, reminiscent of the scrambling time 
for black holes. 

Many further advances are of course necessary before 
one can truly gauge whether the evolution of black holes 
may indeed be understood as the physics of Bose conden¬ 
sates. For one, further improvements of the toy model 
are in line to provide a better understanding of the pro¬ 
cesses in GR. This comprises the inclusion of particle 
number violating vertices and of the entire tower of in¬ 
teractions present in GR, as well as the generalization to 
non-vanishing helicity, and, in hand, the implementation 
of longitudinal modes that are responsible for the grav¬ 
itational potential. Moreover, the compatibility of the 


presence of an instability with the continuously critical 
behavior that seems to be present for black holes needs to 
be understood. This is equivalent to a dynamical under¬ 
standing of the “maximal packing” [5] property of black 
holes, i.e. the fact that it is impossible to localize energy 
beyond its Schwarzschild radius. In the simplest case, 
steps in this direction may be taken already at the level 
of “classicalizing” [26] scalar field theories. It may how¬ 
ever well be that it is the specific structure of General 
Relativity that holds the key to this behavior. 
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Appendix 


In this appendix, we briefly revisit the initial assump¬ 
tions stated in section H. Using our results from section 
VI, we can study their viability by estimating a generic 
decay V —>• N — k. The setup is illustrated in Fig. 6. Due 
to interaction of n constituents, k' particles are emitted, 
while n — k particles reenter the condensate. The grey 
ellipse represents a generic tree-level interaction, yielding 
2n — k + k' — 2 powers of the Planck mass. We require 
that no on-shell fluctuations rescatter on the condensate, 
since such processes receive additional suppression fac¬ 
tors of \/N. Under the assumption that the individual 
momentum transfer scales as (which is the case as 

long as n, k, k' ^ N), we obtain for the squared matrix 



N-n 


FIG. 6: Generic diagram for transition N ^ N — k. n 
constituents interact, leading to the emission of k' 
particles. 
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element 


The squared amplitude is thus bounded by 


\M\^ 


Ar-2K+fc-fc'+2 (-^ ^)- 

{N-n)l {N-n)l 


(A.l) 


which we may approximate by 


\M\^ 


_ 2n+k 

N ) 

/ 2n^ — n + — 

V N 





(A.2) 


IMp < Af2-fc'. (A.3) 

We see that the largest contribution indeed stems from 
the lowest order vertex. Moreover, the scaling from sec¬ 
tion VI survives. 

Note that (A.3) is a crude upper bound on the scaling 
of the higher order interactions. More precise answers, 
in particular on the resummed rate, require an in-depth 
analysis of the specific vertex structure of the theory; in 
the case of GR, additional suppression factors are found 
[ 10 ]. 
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